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THE LOWER BOUNDS FOR EIGENVALUES OF ELLIPTIC OPERATORS 
BY NONCONFORMING FINITE ELEMENT METHODS 



JUN HU*, YUNQING HUANG 1, , AND QUN LIN* 



Abstract. Finding eigenvalues of partial differential operators is important in the mathe- 
matical science. Since the exact eigenvalues are almost impossible, many papers and books 
investigate their bounds from above and below. It is well known that the variational prin- 
ciple (including the conforming finite clement methods) provides the upper bounds, while 
there are no general theories to provide the lower bounds. The aim of our paper is to 
introduce a new systematic method that can produce the lower bounds for eigenvalues. 
The main idea is to use the nonconforming finite element methods. However, the numerics 
from the literature demonstrate that some nonconforming elements produce upper bounds 
of eigenvalues though some other nonconforming elements yield lower bounds. The general 
herein conclusion is that if the local approximation property of the nonconforming finite 
element space Vu is better than the global continuity property of Vh, the corresponding 
method of the eigenvalue problem will produce the lower bound for the eigenvalue. More 
precisely, under three hypothesis on the continuity and approximation properties of the 
nonconforming finite element spaces we first show the abstract error estimates of approxi- 
mate eigenvalues and cigenfunctions. Subsequently, we propose a condition and prove that 
it is sufficient to guarantee the nonconforming finite element methods to produce the lower 
bounds for eigenvalues of the symmetric elliptic operators. As the first application, we 
check the most used nonconforming elements, e.g., the Wilson element, the nonconforming 
linear element by Crouzeix and Raviart, the nonconforming rotated Q\ element by Ran- 
nacher and Turek, and the enriched nonconforming rotated Qi element by Lin, Tobiska 
and Zhou for the second order elliptic operators, the Morlcy element and the Adini element 
for the fourth order elliptic operators, and the Morley-Wang-Xu element for the 2m-th 
order elliptic operator. We show that the condition holds for most of these nonconforming 
elements under some conditions whence they will give lower bounds for eigenvalues of these 
operators. This explains in some sense the numerics in the literature. As the most impor- 
tant application, the condition provides a guidance to modify the known nonconforming 
elements in the literature and to propose new nonconforming elements. In fact, we enrich 
locally the Crouzeix-Raviart element such that the new element satisfies the condition; we 
propose a new nonconforming element for the second order elliptic operators and prove that 
it will yield the lower bounds for eigenvalues. Finally, we prove the saturation condition 
for most of the nonconforming elements aforementioned. As an important application of 
the theory herein, we are able to obtain the upper bound and then design high accuracy 
approximation of the eigenvalue in a very economical way. 
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1. Introduction 

It is a fundamental problem to achieve lower bounds for the eigenvalues of the elliptic 
operators. The reason partially lies in that the Rayleigh-Ritz procedure permits compara- 
tively easy evaluation of upper bounds for such eigenvalues, but the problem of obtaining 
lower bounds is generally considerably more difficult. Moreover, a simple combination of 
the lower and upper bounds will produce an interval to which the exact eigenvalue belongs. 
This in turn gives a reliable a posteriori error estimate of the approximate eigenvalue, which 
is essential for the design of the coefficient of safety in practical engineering. In fact, the 
study of the lower bound for the eigenvalue can date back to the remarkable works of [331 El] 
and [88, 89J, where the lower bounds of the eigenvalues are derived by the finite difference 
method for the second order elliptic eigenvalue problem. The method based on the in- 
termediate problems [9j [TlJ [101 [121 EDI EQl [35] gives upper and lower bounds for certain 
domains. We refer the interested reader to the paper [35] for a comprehensive survey on 
the method of intermediate problems in the earlier time. Moreover, the method of collo- 
cation by interpolation in combination with the method of particular solutions admits the 
approximate eigenvalue from below. The homotopy method can also produce lower bounds 
for eigenvalues [HI H2j [72]. If we only consider the first eigenvalue, we can refer to a very 
wonderful method proposed in [70]. We also refer the interested readers to [54J for the 
various numerical methods for the eigenvalues of the Laplacian operator in two dimensions. 
However, the aforementioned constructive methods depend on both the operators and the 
domains as well as the boundary conditions. In particular, they will become complicated 
for the complex domain. 

Since that the finite difference method in some sense coincides with the standard linear 
finite element method with mass lumping, one could expect that the finite element method 
with mass lumping gives lower bounds for eigenvalues of the operators. For the uniform 
mesh, such a result is actually proved in [16| Wf\ . For the general case, it is shown that the 
approximate eigenvalues produced by the linear finite element method with mass lumping 
are always smaller than those by the standard linear finite element method [IHl H7J, see 
also [3] . However, the numerical examples in both [161 Wf\ and [3] suggest that the mass 
lumping method admit the lower bounds for the eigenvalues provided that the eigenfunction 
is smooth, see [60] for the proof of the one dimensional case. 

The nonconforming finite element methods are alternative possible ways to produce lower 
bounds for eigenvalues of the operators. In fact, from the very beginning, engineers believe 
that the incompatibility of the nonconforming methods makes the structure less stiff, see 
for instance, [S. M. Lepi, Practical Guide to Finite Elements: A solid mechanics approach, 
Marcel Dekker, Inc, Page 209, 1998; A. Ghali, A. M. Neville and T. G. Brown, Structural 
analysis: a unified classical and matrix approach, Spon Press, London, Page 526, 2003] 
and the references therein. Mathematicians also think that both the strain energy and the 
displacement converge from above. For example, Strang and Fix remark for the noncon- 
forming Wilson element that " We are confident that this is the correct rate of convergence, 
although we must emphasize two points: (1) It is obvious experimentally that the constant 
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which multiplies h is very much smaller than it would be without the extra nonconforming 
trial functions, and (2) the strain energy in is no longer necessarily smaller than that in 
u (where is the nonconforming finite element approximation to u and h is the meshsize of 
the mesh). In fact, convergence from above is the rule rather than the exception, both for 
strain energy and for displacements", see, [G. Strang and G. Fix, An Analysis of the Finite 
Element Method, Prentice-Hall, Page 180, 1973]. Although there are no theoretical justifi- 
cation for the remarks from both engineers and mathematicians, the lower bound properties 
of eigenvalues by nonconforming elements are observed in numerics, see, Zienkiewicz et al. 
[lOlj . for the nonconforming Morley element, Rannacher [73], for the nonconforming Mor- 
ley and Adini elements, Liu and Yan [68], for the nonconforming Wilson [931 EH E2], EQ™* 
[64] . and Qi ot UM elements. See, Boffi [16], for further remarks on the possible properties 
of discrete eigenvalues produced by the nonconforming methods. 

There are a few results to study the lower bound property of eigenvalues by nonconforming 
elements. The first result in this direction is analyzed in a remarkable paper Armentano and 
Duran [I] for the Laplacian operator without a lower order term on the polygonal domain 
when the homogeneous Dirichlet boundary condition is imposed. The analysis is based on 
an identity for the error of the eigenvalue. It is proved that the nonconforming linear finite 
element method of [25] leads to the lower bounds for eigenvalues provided that the eigen- 
function u G H 1+r (Q) (1 Hq(Q) with < r < 1. In some sense, this means the approximate 
eigenvalue will be smaller than the exact one if its associated eigenfunction is singular. The 
idea is generalized to the Wilson element in Zhang, Yang and Chen [100J, where the identity 
for the error of the eigenvalue is first developed to the general nonconforming finite element 
methods. At the same time, the method of [I] is extended to the enriched nonconforming 
rotated Qi element of [M] in Li [58J. It is proved that the approximate eigenvalues are 
always smaller than the exact ones provided that u G H 1+r (Q) R Hq(Q) with < r < 1. 
Note that the case with r = 1 is allowed as well therein. Therefore this covers both the case 
with singular eigenfunctions and the case with smooth eigenf unctions. A similar result is 
analyzed for the enriched nonconforming rotated Qi element in [97] where the nonconform- 
ing linear element is also investigated. The extension to the Morley element can be found 
in [HE]. However, all of those papers except |100j depend on the saturation condition for 
which a rigorous proof is missed in the literature. 

The expansion method based on superconvergence or extrapolation is also used to analyze 
the lower bound property of eigenvalues by the nonconforming elements [100l[97 ][65l[63| [62]. 
The research in this direction can date back to a series of talks "Finite Element Method: 
Accuracy and Extrapolation" by the third author in Hebei University in 2001. In these 
talks, he derived the expansions and extrapolations of errors of the eigenvalues by the 
nonconforming rotated Qi element, the Wilson element and the enriched nonconforming 
rotated Q\ element on the uniform rectangular meshes. These results were later written into 
the book [65] . The numerical tests of these results can be found in [671 EE] . The lower bound 
property of eigenvalues by the Wilson element in two dimensions is also analyzed in [631 1100[ 
[97] , while the lower bound property of eigenvalues by the enriched nonconforming rotated 
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Qi element is studied in [621 [9?]. It is proved in [951 EH] that the approximate eigenvalues 
produced by the nonconforming Adini element are smaller than the exact solutions provided 
that the meshsize is small enough. Most of the aforementioned results are summarized in a 
recent survey paper [97]. We note that the expansion method is only valid for the smooth 
enough case on the uniform rectangular mesh. 

The aim of our paper is to introduce a new systematic method that can produce the lower 
bounds for eigenvalues. The main idea is to use the nonconforming finite element methods. 
However, the numerics from the literature demonstrate that some nonconforming elements 
produce upper bounds of eigenvalues though some other nonconforming elements yield lower 
bounds, see[73, 68J. We find that the general condition lies in that the local approximation 
property of the nonconforming finite element space Vh should be better than the global 
continuity property of V^. Then the corresponding nonconforming methods will produce 
the lower bounds for eigenvalues of the elliptic operators. More precisely, first, we shall 
analyze the errors of the discrete eigenvalues and eigenf unctions. We refer the interested 
readers to [15l ESI HZ1 EH EH EHl EHl E21 EH [71 El [521 H3] for the convergence analysis in 
other setting. Second, we shall propose a condition on the nonconforming element methods 
and then under the saturation condition prove that it is sufficient for the lower bounds for 
eigenvalues. With this result, to obtain the lower bound for the eigenvalue is to design 
nonconforming finite element spaces with enough local degrees of freedom when compared 
to the global continuity. This in fact results in a systematic method for the lower bounds of 
eigenvalues. As the first application of our method, we check that this condition holds for the 
most used nonconforming element, e.g., the Wilson element [93 1 1771 179] , the nonconforming 
linear element by Crouzeix and Raviart [25j, the nonconforming rotated Qi element by 
Rannacher and Turek [7H [79] , and the enriched nonconforming rotated Q\ element by Lin, 
Tobiska and Zhou [64J for the second order elliptic operators, the Morley element [69l [79] 
and the Adini element [79] for the fourth order elliptic operators, and the Morley- Wang- 
Xu element [H7] for the 2m-th order elliptic operator. These recover all the results in 
[U EH EH1 11001 1971 1651 1631 162] . As the most important application, we follow this guidance 
to enrich locally the Crouzeix-Raviart element such that the resulted element satisfies the 
sufficient condition and to propose a new nonconforming element method for the second 
order elliptic operators and show that it actually produces the lower bounds for eigenvalues. 
As the indispensable and important part of the paper, we prove the saturation condition 
for most of these nonconforming elements under consideration. Compared to the results in 
the literature, the new contributions of the paper are as follows: 

• Above all, we introduce a new systematic method that can produce lower bounds 
for eigenvalues of symmetric elliptic operators. The users can follow the condition 
and examples herein to construct new nonconforming methods that can produce the 
lower bounds for eigenvalues of operators under their consideration. The advantage 
of the method over other methods in the literature is that the method herein is 
independent of the operators and the boundary conditions. Moreover, these lower 
bounds are available up to the generalized eigenvalue problems of matrixes. 
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• We analyze the nonconforming Wilson element method for the general second order 
elliptic eigenvalue problems with general boundary conditions in any dimension on 
nonuniform n-rectangular meshes and prove that the discrete eigenvalues are smaller 
than the exact ones for both the singular and smooth cases. Moreover, as we point 
out below that it is straightforward to extend the analysis and result to the general 
quadrilateral mesh in two dimensions and the general hexahedron mesh in the three 
dimensions. 

• We follow the sufficient condition to propose two classes of nonconforming methods 
that produce the lower bounds for eigenvalues of the second order elliptic operators. 

• We analyze the nonconforming Adini element method for the general fourth order 
elliptic eigenvalue problems with the general boundary conditions in any dimension 
on nonuniform meshes and prove that the discrete eigenvalues are smaller than the 
exact ones for the smooth case. 

• We analyze the Morley-Wang-Xu element for the 2m-th order elliptic eigenvalue 
problems and prove that the discrete eigenvalues are smaller than the exact ones for 
both the singular and smooth cases. 

• We propose and prove a sufficient condition for the saturation condition of the piece- 
wise polynomial approximation and consequently establish the saturation condition 
for most of nonconforming elements under consideration. 

As an important application of the theory herein, we are able to obtain the upper bound 
and then design high accuracy approximation of the eigenvalue in a very economical way. 
In fact, we can get a conforming approximation by some simple averaging interpolation [18] 
of the nonconforming approximation of the eigenfunction. Now the Rayleigh quotient of 
the resulted conforming approximation will be the upper bound of the eigenvalue under 
consideration. Then we can combine the lower bound and the upper bound to achieve high 
accuracy approximation of the eigenvalue [49] . 

The paper is organized as follows. In the following section, we shall present the symmet- 
ric elliptic eigenvalue problems and their nonconforming element methods in an abstract 
setting. In Section 3, based on three conditions on the discrete spaces, we analyze the error 
estimates for both the discrete eigenvalues and eigenfunctions. In Section 4, under one more 
condition, we prove an abstract result that the eigenvalues produced by the nonconforming 
methods are smaller than the exact ones. In Sections 5-7, we check these conditions for 
various nonconforming methods in the literature and we also propose two new classes of 
nonconforming methods that admit the lower bounds for eigenvalues in Section 8. We end 
this paper by Section 9 where we give the conclusion, which is followed by the appendix 
where we analyze the saturation condition for the piecewise polynomial approximation. 

2. The eigenvalue problems and nonconforming finite element methods 

Let V C H m (Q) denote some standard Sobolev space on some bounded Lipschitz domain 
fl in M n with a piecewise flat boundary dfl. The 2m-th order elliptic eigenvalue problem 
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reads: Find (X,u) G M x V such that 

(2.1) a(u,v) = \{pu,v) L 2( Q) for any v G V and \\p 1/2 u\\ L 2 {n) = 1, 

with some positive function p G L°°($7). The bilinear form a(u,v) is symmetric, bounded, 
and coercive in the following sense: 

(2.2) a(u, v) = a(v, u), \a(u, v)\ < ||rx||y||u ||y, and \\v\\ v < a(t>, t>) for any u,v E V, 

with the norm \\-\\y over the space 1/. Throughout the paper, an inequality A < B replaces 
A < C B with some multiplicative mesh-size independent constant C > that depends only 
on the domain Q, the shape (e.g., through the aspect ratio) of elements, and possible some 
norm of eigenfunctions u. Finally, A « B abbreviates A < B < A. 

Under the conditions (12. 2p . we have that the eigenvalue problem (12.11) has a sequence of 
eigenvalues 

< Ai < A 2 < A 3 < ■ • • / +oo, 
and corresponding eigenfunctions 

(2.3) ui,u 2 ,u 3 , ■■■ , 
which can be chosen to satisfy 

(2.4) (pUi,Uj) L 2 {n) = 5 i:j ,i,j = 1,2, • • ■ . 
We define 

(2.5) E t = span{ui,« 2) • • • ,ui}. 

Then, the eigenvalues and eigenfunctions satisfy the following well-known minimum-maximum 
principle: 

. a(v,v) a(u,u) 
(2.0) \k = mm max r = max ■ 



For any eigenvalue A of (12.11) . we define 

(2.7) M(A) := \u : u is an eigenfunction of (12.11) to A}. 

We shall be interested in approximating the eigenvalue problem (12.11) by finite element 
methods. To this end, we suppose we are given a discrete space Vh defined over a regular 
triangulation Th of Q into (closed) simplexes or n- rectangles [19]. 

We need the piecewise counterparts of differential operators with respect to Th- For any 
differential operator £, we define its piecewise counterpart Ch in the following way: we 
suppose that vk is defined over K 6 Th and that the differential action Cvk is well-defined 
on K which is denoted by CrVr for any K e %] then we define Vh by Vh\x — V K where 
Vh\x denotes its restriction of Vh over K\ finally we define ChVh by {ChVh)\K = C-kVr- 

We consider the discrete eigenvalue problem: Find (A/j, Uh) Gf xV ft such that 

(2.8) a h (u h ,v h ) = \h{pUh,v h ) L 2(n) for any v h G V h and \\p 1/2 u h \\ L 2 (n) = 1. 
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Here and throughout of this paper, a/ l (-, •) is the piecewise counterpart of the bilinear form 
a(-, •) where the differential operators are replaced by their discrete counterparts. We assume 
that the discrete space Vh satisfies 

(HI) || • := ah(-, -) l l 2 is a norm over the discrete space Vh- 

Let N = dimVft. Under this condition, the discrete problem (12. 8 1) admits a sequence of 
discrete eigenvalues 

< Xx t h < A2,/i < • • • < Ajv,ft, 
and corresponding eigenfunctions 

V>1, h, U2,h, • • " j u N,h ■ 

In the case where Vh is an conforming approximation in the sense Vh C V, it immediately 
follows from the minimum-maximum principle (12. 6 j) that 

Afe < A*.,/!, fc = 1,2, • • ■ , iV, 

which indicate that A^ is an approximation above to A&. 
We define the discrete counterpart of Eg by 

(2.9) E tjh = span{«i ih , u 2 , h , • • • , u^ fc }. 
Then, we have the following discrete minimum-maximum principle: 

/ 01f y\ \ • a h (u,u) 

(2.10) Afc^ = mm max = max 



dimV k<h =k,V k<h cV h veV hlh (pV,v) L 2^ u£E k}h (p«, w)l 2 (Q) 

3. Error estimates of the eigenvalues and eigenfunctions 

In this section, we shall analyze the errors of discrete eigenvalues and eigenfunctions 
by the nonconforming methods. We refer to [521 [75] for some alternative analysis in the 
functional analysis setting. We would like to stress the analysis is a nontrivial extension 
to the nonconforming methods of the analysis for the conforming methods in [81] . For 
simplicity of presentation, we only consider the case where Xi is an eigenvalue of multiplicity 
1. 

Associated with the bilinear form a(-, •), we define the operator A by 

(3.1) a(u, v) = (Au,v) L 2rm for any v G V . 

The conditions on the approximation and continuity properties of the discrete space Vh 
are assumed as follows, respectively, 

(H2) Suppose v G V D H 3+ (ft) with some integer j > m and < S < 1. Then, 

inf \\v - v h \\ h < h 3+s - m \v\ w+ s m . 

(H3) Suppose v G V D H h+S (Q) with some integer k > m and < s < 1. Then, 

a h (v,Vh) - (Av,v h ) L 2 {Q) 
sup 7, — fi ^ ft |w|H fe + s (n) , 

with m <k < ] and s < 5 when k — j. 
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Given any / G L 2 (Q), let Uf be the solution to the dual problem: Find Uf G V such that 
(3.2) a(u f , v) = (pf, v) L 2 {n) for any v G V . 

Generally speaking, the regularity of Uf depends on, among others, the regularities of / 
and p, the elliptic operator under consideration, the shape of the domain Q and the bound- 
ary condition imposed. To fix the main idea and therefore avoid too technical notation, 
throughout this paper, without loss of generality, we assume that u f G V H H k+S (tt) in the 
sense that 

(3-3) \\u f \\ H k +sm < \\p 1/2 f\\ L 2(n) ■ 

In order to analyze the L 2 error estimate of the eigenfunction, we define the quasi-Ritz- 
projection P' h ug G V h by 

(3.4) ah{P'h u ti v h) = >>e(pug, v h ) L 2 (n) for any v h G V h . 

Throughout this paper, ug, Uj, and Ui are the eigenfunctions of the problem (12.11) . while 
u e,h, Uj,h, an d u i,h are the discrete eigenfunctions of the discrete eigenvalue problem. Note 
that P' h ug is the finite element approximation of ug. Under conditions (H1)-(H3), a standard 
argument for the nonconforming finite element method, see, for instance, [2U[T9l[79], proves 

Lemma 3.1. Suppose that the hypothesis (Hl)-(HS) hold. Then, 

(3.5) Wp^iug - PLug)\\ L 2 {n) + h k+s ~ m \\ug - P' h ug\\ h < h^ k+s -^\ug\ Hk+ s {n) . 
From P' h ug G Vh we have 

N 

(3.6) P' h ug = ^{pP' h Ug : Uj th )u jjh . 

For the projection operator P h , we have the following important property 

(3.7) {\j,h - Xg){pP' h Ug,U hh ) L 2 {n) = \g(p(ug - P' h Ug),Uj. h ) L 2 {n) . 

In fact, we have 

(3.8) \ jt h{pP' h ug, Uj,h)L*(ci) = ah(uj,h, PhUg) = Xg(pug, Uj,h)&(ci) ■ 
Suppose that ^ Xj if t ^ j. Then there exists a separation constant dg with 

(3.9) Xl < d e for any j ^ £, 
provided that the meshsize h is small enough. 

Theorem 3.2. Letug andug^ be the eigenfunctions of (12. ip and (12. 8p . respectively. Suppose 
that (HI)- (H3) hold. Then, 

(3.10) \\p l ' 2 (ug - u e , h )\\ L 2 {a) < h 2( ~ k+s -^\ug\ Hk+s{n) . 
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Proof. We denote the key coefficient {pP' h ug^ue,^) by fa. The rest can be bounded as 
follows: 

\\p 1/2 (P^U e - fau ith )\\ 2 L 2 in) = ^2(pPLu e ,U jA ) 2 L 2 (n) < 4 ^2(P(UI - PLut),U jih ) 2 L2{n) 

(3.11) m a* 

<4\\p Xl \^-P' h Ul)f L ^y 

This leads to 

\\p 1/2 (u e - faui th )\\ L 2 {n) < \\p 1/2 (u e - PhUi)\\i?(fl) + \\p 1/2 (Pl l u i - faui )h )\\ L 2( a) 
(3-12) , w 

< (1 + d^iut - PMWvw $ h 2 ^ s - m hf +s)/2m . 

\\p 1/2 ue\\L2(n) ~ \\p 1/2 {u£ - fau eth )\\ L 2 {n) < \\p 1/2 fau eih \\ L 2 (n) 

(3.13) 

< \\p 1/2 u e \\ L 2( n) + \\p 1/2 (ue - faut h )\\ L 2 {n) . 
Since both ug and u^h are unit vectors, we can choose them such that fa > 0. Whence we 
have \fa - 1| < \\p 1/2 {u e - IU a (n)- Thus, we obtain 

\\p l,2 {ut - uth)\\msi) < \\p 1/2 (ui - fau ijh )\\ L 2( U) + \fa - l|||p 1/2 M^||i,2m) 

(3.14) 

< 2\\p 1 / 2 (u e - fau^ h )\\ L 2 m < h 2 ( k+s ^\u e \ Hk+s(n) . 
This completes the proof. □ 

Next we analyze the error of the eigenvalue. To this end, we define u^h £ V by 
(3.15) a(u eih ,v) = \h{pUL h ,v) L 2 [n) for any v E V. 

It follows from (J2J) and (EHD that 

(p(ue,h - ut,h),ui)L 2 (Q) = Xt l \h{put,h,ut) L *(ti) - (pUth,Ui) L 2( Q) 

Thus we have 

» » Mp(ui,h-u e , h ),u e ) L 2 {n) 
[o.it ) M,h~M — 7 s • 

{pU £)h ,U £ ) L 2 {n) 

For simplicity of notation, we assume that u^h G V" fl H k+S (Q) with < s < 1 in the sense 
that 

(3-18) ^ ^e,h\\p 1 ^ 2 ue,h\\L 2 (n)- 

Note that u^h is the finite element approximation of u^h. A standard argument for the 
nonconforming finite element method, see, for instance, [2"l"l [T9"| [75]. proves 

Lemma 3.3. Suppose that the hypothesis (Hl)-(HS) hold. Then, 

(3.19) \\p 1/2 (u e , h - UL h )\\ L 2 (n) + h k+s - m \\u eA - u eA \\ h < X e , h k+s -^\\p^ 2 u e>h \\ L 2 {n) . 

Inserting the above estimate into ( 13.171) proves: 
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Theorem 3.4. Let \e and A^ be the eigenvalues of (12.11) and (12.81) . respectively. Suppose 
that (H1)-(H3) hold. Then, 

(3.20) IV - A/| < h 2 ^ k+s -^\u e \ Hk+ s {n) , 
provided that h is small enough. 

Finally we can have the error estimate in the energy norm of the eigenfunction. 

Theorem 3.5. Letui andu^h be the eigenf unctions of (12. ip and (12. 8p . respectively. Suppose 
that (HI)- (H3) hold. Then, 

(3.21) \\u t - u iih \\ h < h k+s - m \u e \ Hk+ s {n) . 

Proof. In order to bound the error of the eigenfunction in the energy norm, we need the 
following decomposition: 

a h (u e - u Ql , u £ - ue >h ) = a(ut, u e ) + a h (u £th , u £jh ) - 2a h (u e , u e>h ) 

(3.22) 

= \i\\p 1/2 (ue - ue,h) \\h(n) + \h - A £ + 2Xi(pui, u i<h - u e ) - 2a h {u£, u^ h - ug) . 
Then, the desired result follows from Theorem 13.41 (I3.14p . and the hypothesis (H3). □ 
We end this section to define the Galerkin projection operator : V — > 14 by 

(3.23) a h (P h v, w h ) = a h (v, w h ) for any w h e V h , v e V. 

Remark 3.6. We note that P' h will become Ph for the conforming methods, which indicates 
the difference between the conforming elements analyzed in [HI] and the nonconforming 
elements under consideration. 

Under conditions (HI), (H2), and (H3), a standard argument for the nonconforming finite 
element methods, see, for instance, [2U [HI ES] , proves 

(3.24) ^(v - P h v)\\ LHn) + h k+s - m \\v - P h v\\ h < h 2 ^ k+s -^\v\ Hk+ s {n) , 
provided that v e V n H k+S {Q). 

4. Lower bounds for eigenvalues: an abstract theory 

This section proposes a condition on the nonconforming element method and proves that 
it is sufficient to guarantee the nonconforming method to yield lower bounds for eigenvalues 
of the operators. 

In the following, we shall use some identity for the error of the eigenvalue. This idea is 
first proposed in Armentano and Duran[l], which is developed in Zhang, Yang, Chen [100J, 
Li [5E], Yang, Lin, and Zhang [96], Yang, Zhang, and Lin [97], and Yang, Lin, Bi, and Li 
[98] . The identity with the full terms first appears in Zhang, Yang and Chen [100] . We shall 
give the full details for readers' convenience. 
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Lemma 4.1. Let (X,u) and (\h,Uh) be the solutions of problems ( 12. lft and ( 12. 8p . respec- 
tively. For any Vh G V^, we /lave t/ie following identity: 

X-X h = \\u - u h \\ 2 h - X h \\p 1/2 (v h - u h )\\ 2 L2{n) 

(4.1) 

+ A h (||p 1/2 u h ||ia ( n) - \\p 1/2 u\\h(n)) + 2a h (u-v hl u h ). 

Proof. From the eigenvalue problem (12.11) and the discrete eigenvalue problem (12.81) we have 

X + X h = \\u\\l + \\u h \\l 

= \\u-u h \\l + 2a h (u,u h ) 

= \\u - u h \\ 2 h + 2a h (v h , u h ) + 2a h (u - v h , u h ) 

= \\ u ~ u h\\l + 2\ h {pv h ,u h ) L 2 i p j) + 2a h (u - v h ,u h ). 

This leads to 

A - X h = \\u - Uh\\l + 2X h (p(v h - u h ),u h ) L 2( Q) + 2a h (u - v h , u h ) 

= \\u - U h \\ 2 h - X h (p(v h - U h ),V h - U h ) L ?(Q) 

+ X h (p(v h - Uh), v h + Uh)i?(fi) + 2a h (u - v h , u h ) 

= \\u - U h \\l - X h (p(v h - U h ), V h - Uh)L2(Q) 

+ Xh{\\p 1/2 v h \\ 2 L 2 {n) - \\p 1/2 u\\ 2 L 2 {n) ) +2a h (u-v h ,u h ), 
which completes the proof. □ 

The essential requirement on the discrete space Vh is hidden in the following condition: 
(H4) Let u and Uh be the eigenfunctions of problems (12.11) and (12.8p . respectively. We 
assume that there exists an interpolation TlhU G Vh with the following properties: 

a h {u - U h u,u h ) = 0, 
(4.2) |llp 1/2 -lli 2( ,) - \\p l/2 n h u\\ 2 L 2 (n) \ < h***™***, 

when the meshsize h is small enough and u G V fl H J+ (Q) with m < k < j, < «S < 1, 
< s < 1 and s < S when k — j, two constants < As and < AS. 
From the abstract error estimate (I3.10p we have 

(4-3) \\p 1/2 (v-Uh)\\vm <h 2{k+s ~ m K 

Whence the triangle inequality and (H4) shows that the second and third terms on the 
right-hand side of (14.11) are of higher order than the first term. Since = ah{u — UhU,Uh), 
we have the following main result: 

Theorem 4.2. Let (X,u) and (Xh,Uh) be the solutions of problems (12.11) and (I2.8p . re- 
spectively. Assume that u G Vfl G H' J+S {VL) with m < k < j and < S < 1 and that 
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^2{k+s m) < || n _ y^^/j < s < 1 and s < S when k = j. If the four assumptions (HI) 
-(H4) hold, then 

(4.4) \ h < A, 

provided that h is small enough. 

The condition that h 2 ( k + s - m ) < \\u — Uh\\h is usually referred to as the saturation condition 
in the literature. The condition is closely related to the inverse theorem in the context of the 
approximation theory by the trigonometric polynomials or splines. For the approximation 
by the conforming piecewise polynomials, the inverse theorem was analyzed in [921 E| • For 
the nonconforming finite element methods, the saturation condition was first analyzed in 
Shi |78j for the Wilson finite element method through an example, which was developed by 
Chen and Li [23] by an expansion of the error. See [99] for some example of the conforming 
linear element in the two dimensions, see [53] for the lower bound of the discretization error 
by the conforming linear /bilinear finite element. Babuska and Strouboulis [6] analyzed 
the Lagrange finite element method for the elliptic problem in the one dimension. In the 
appendix, we shall analyze the saturation condition for most of nonconforming finite element 
methods under consideration. To our knowledge, it is the first time to analyze systematically 
this condition for the nonconforming methods. In particular, we first prove the saturation 
condition when the eigenfunction is singular in the sense u G H m+S (Q) with < s < 1 by 
combining the interpolation space theory and the idea from a posteriori error analysis of 
the nonconforming finite element methods developed recently by Hu and Xu |50j . 

We will take the Galerkin projection operator Ph from (13. 23ft or its high order perturbation 
not the canonical interpolation operator of the nonconforming space Vh as the interpolation 
operator Uh- Then the error estimate of the operator Ph is dependent on only the local 
approximation property but not the global continuity property of the space Vh while the 
error — Uh\\h generally depends on both properties (see Theorems 13.21 and 13.51 ). On the 
other hand, the term Hp 1 ^ 2 ^!!^^) — Hp^n/i^H 2 ^^) will be of high order when we have enough 
many local degrees of freedom (compared to the continuity) and therefore the consistency 
error in \\u — W/i|| 2 will be dominant in the sense that 

\\u-u h \\l > Il|p 1/2 w|li2 (n) - ||p 1/2 n h u||£a (n) |. 

If this happens we say the local approximation property of the space Vh is better than 
the global continuity property of Vh- Whence, Theorem 14.21 states that the correspond- 
ing method of the eigenvalue problem will produce the lower bound for the eigenvalue for 
this situation. Thus, to get a lower bound for the eigenvalue is to design nonconforming 
finite element spaces with enough local degrees of freedom when compared to the global 
continuity property of Vh- This in fact provides a systematic tool for the construction of 
the lower bound for the eigenvalue of operators in mathematical science. We would like to 
stress that the theory herein applies to all symmetric elliptic operators not restricted to the 
2m-th symmetric elliptic operators under consideration; cf. [66] for the extension to the 
nonconforming elements for the Stokes problem. 
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In the subsequent sections, we shall first show that the condition (H4) expressed in ( 14.2ft 
holds for most of nonconforming finite element methods in the literature. Then we shall 
modify some nonconforming elements which does not satisfy the condition (H4). We prove 
that the modified method will meet the condition (H4) and consequently they will produce 
the lower bound of the eigenvalue. Finally, we follow the condition (H4) to propose some 
new nonconforming finite element method which yields the lower bound of the eigenvalue. 

5. Nonconforming elements of the second order operators 

In this section, we shall present some nonconforming schemes of the second order elliptic 
eigenvalue problem that the conditions (H1)-(H4) proposed in the previous sections are 
satisfied. Let the boundary dfl be divided into two parts: and with |ro| > 0, and 

r D u r N = an, and r D n r N = 0. 

The eigenvalue problem reads: Find (A, u) G R x V =: R x H]j{Q) —: R x {v G 
H\tt),v\r D = 0} with 

(5.1) a(u,v) L 2 {n) = X(pu,v) L 2 {Q) for any v G V with \\p 1/2 u\\ L 2 {n) = 1, 

where p G L°°(Q) is a positive function. In this section we only assume that (15.11) is a 
general second order elliptic eigenvalue problem imposed the general boundary conditions 
unless otherwise stated. In particular, we assume that the bilinear form a(-, •) satisfies the 
conditions of (12.21) . 

Let Th be a regular n-rectangular triangulation of the domain Q C R n with 2 < n in 
the sense that U^e^ ^ = ' ^ wo distinct elements K and K' in Th are either disjoint, or 
share the ^-dimensional hyper-plane, t = 0, • • • ,n — 1. Let Hh denote the set of all n — 1 
dimensional hyper-planes in Th with the set of interior n — 1 dimensional hyper-planes %h{&) 
and the set of boundary n — 1 dimensional hyper-planes 1-LhipVi). Afh is the set of nodes of 
Th with the set of internal nodes and the set of boundary nodes Mh(dQ). 

For each K G Th, we introduce the following affine invertible transformation 

F K : K -> K,Xi = h Xi)K £i + x°i 

with the center (x?, x®, ■ • ■ , x°) and the lengths 2h Xi) K of K in the directions of the x^-axis, 
and the reference element K = [—1, l] n . In addition, we set h = maxi<j<„ h Xi . 

Over the above mesh Th, we shall consider three classes of nonconforming element methods 
for the eigenvalue problem (15.11) . namely, the Wilson element in any dimension, the enriched 
nonconforming rotated Q\ element in any dimension, and the nonconforming rotated Q\ 
element in any dimension in the literature. 

Let Vh be the discrete spaces of the aforementioned nonconforming element methods. The 
finite element approximation of Problem (15.11) reads: Find (Xh,Uh) G R x Vh, such that 

(5.2) a h (u h ,v h ) L 2( Q) = \h{pu h ,v h ) L 2(n) for any v h G V h and \\p 1/2 u h \\ L 2 {n) = 1 , 

with the piecewise defined counterpart ah(-, •) of a(-, •). 

For all the elements, one can use the continuity and boundary conditions for the discrete 
spaces V h given below to verify the hypothesis (H1)-(H3), see [931 El E3 [79] for 
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further details. Let (A, it) and (A^, Uh) be the solutions to the problems (15. ip and (15. 2D . by 
Theorems 13.21 13. 4[ and 13.51 we have 

(5.3) |A - X h \ + h s \\u - u h \\ h + \\p 1/2 (u - u h )\\ L 2 {n) < h 2s , 

provided that u G V n # 1+0 (fi) with < s < I. 

We shall analyze the key condition (H4) for these elements in the subsequent subsections. 

5.1. The Wilson element in any dimension. Denote by Qwu{K) the nonconforming 
Wilson element space [931 EH EH] on the reference element defined by 

(5.4) Qwu(K) = Q^K) + span{£ 2 - 1, £ 2 2 - 1, • • • , £ ~ 1}, 

where Q±(K) is the space of polynomials of degree< 1 in each variable. The nonconforming 
Wilson element space Vh is then defined as 

Vh '■= { v £ L 2 (Q) : v\k ° Fk £ Qwu{k) for each K G Th,v is continuous 
at the internal nodes, and vanishes at the boundary nodes on To} . 

In order to show the hypothesis (H4), we let be the Galerkin projection operator defined 
in (13.231) . Then, we have 

(5.5) \\u-P h u\\ h < inf \\ u -v h \\ h <h l+ * , 

Vh&Vh 

provided that u G V fl H 2+S (Q) with < s < 1. In order to derive the L 2 error estimate of 
u — PhU, we introduce the following auxiliary problem: Find w G V such that 

(5.6) a(w,v) = (p(u - Ph,u),v) L 2(ti) for any v G V . 
We continue the assumption in (13. 2p to assume that 

(5-7) |M|ip(n) ~ \\P 1/2 ( U ~ Phu)\\ L 2(n) , 

for the case u G V fl H 2+S (Q). Then we have 
(p(u - P h u), u - P h u) L 2i n) 

(5.8) 

= (p(u - P h u),u - P h u) L 2 (n) - a h (w, u - P h u) + a h (w - P h w, u - P h u). 

The first two terms on the right-hand side of the above equation are the consistency error 
term for the Wilson element. Then, we have 

(5-9) \\p l l 2 {u-P h u)\\ L 2 {n) <h 2+s . 

To analyze the singular case where u G V fl H 1+S (Q) with < s < 1, we need the following 
conforming space 

Vh := { v e H 1 ^) '■ v\k ° Fk G Qi(K) for each K G 7h, v is continuous 

at the internal nodes, and vanishes at the boundary nodes on ^d} ■ 

Given Vh G Vh, we define the interpolation H c h Vh G V£ by 

(5.10) (U c h v h ){V) = v h (V) for any vertex V of % ■ 
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Given v G L 2 (K), we define J' K v G [w o F K G span{^ — 1, • • • , ^ — 1} j such that 

(5.11) / J' K vdx = I vdx. 

Jk Jk 

Given any v G V, we define Jrv G Qwu(K) by 

(5.12) J K v\ K = (U c h P h v)\ K + J' K (v - (U c h P h v)\ K ) . 
The interpolation J^v G Vh is defined by 

(5.13) J h v\ K = J k v\k for any K G 7^. 

For the interpolation operator J^, we have the following crucial conservative property 



(5.14) / Jhvdx = \ vdx for any K e %■ 

Jk Jk 

Lemma 5.1. For the Wilson element defined in (15.51) . there holds the hypothesis (H4) with 
\\p 1/2 (u - P h u)\\ L 2 (n) < h 2+s ifu G H^(Q) n H 2+S (n) with < 5 < 1 , 

(5.15) 

Wp 1 ' 2 ^ - P h u)\\ L 2 {n) < h 1+s ifue Hhfp.) n H 1+S (ty with o < s < l . 

Proof. Let LT^ = P ft in (H4). Since a h {u — P^, Vh) = for any Vh G Vh, we only need to show 
the two estimates. 

For the case where u G H\)(Vt) fl H 2+S (Q) with < s < 1, the estimate (15. 9ft proves two 
estimates in (H4) by k — j — 2, m — 1, s = 0, S — s, As = s, and A5 = 1. 

To prove the result for the case where u G H]j(Q) HH 1+S (Q) with < s < 1, we shall use 
the interpolation operator Jh- Given K G Tk, JhU — PhU vanishes at all the vertexes of K. 
By the scaling argument and the equivalency of the norm of the finite element space on the 
reference element K, this yields 

(5.16) \\J h u - P h u\\ L 2 {K) < h K \\V(J h u - P h u)\\ L 2 {K) . 

By virtue of the local conservative property of J^u in (I5.14p . we use the Poincare inequality 
to get 

||W - PhU\\ L 2(K) < \\J~hU - PhU\\ L 2(K) + \\u - JhU\\tf((l) 

(5.17) 

< h K \\V(J h u - Phu)\\ L 2(K) + h K \\V{u - J h u)\\ L 2 {K ) . 

By the fact that JhV = v for any v G V h c , and the interpolation space theory [T9(, HI] , we 
have 

(5.18) ||u - Phu\\v(n) < h 1+s \u\ H i+ s(n) when u G H l D {Vl) n H 1+5 {Vt) with < s < 1 , 

which proves the two estimates in (H4) by k = m = j = 1, s = S, As = 1 — s, S = s, and 
AS = 1. □ 

In the appendixes B and C , we prove that h < \\u — Uh\\h when u G H 2+S (Q) with 
< s < 1 and that there exist meshes such that h s < \\u — Uh\\h holds when u G H 1+B (Q) 
with < 5 < 1. In view of the estimate ( 15.31) . it follows from the interpolation error estimate 
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(15.15!) that the result in Theorem 14.21 holds for the Wilson element method of the eigenvalue 
problem in any dimension. 

For the Laplace operator, the Wilson element method of the eigenvalue problem in two 
dimensions with n = 2 has already been analyzed in [63]. 11001 [97] under the regularity 
condition u G (Q) R H 3 (Q); the Wilson element method of the eigenvalue problem in any 
dimension has already been analyzed in [61J under the regularity condition u G Hq(Q) fl 
H 2+S (Q) with < s < 1. The analysis of these papers is based on the canonical interpolation 
operator of the Wilson element. However, these results heavily depend on the uniformity of 
the mesh and the homogeneous Dirichlet boundary condition. Therefore, the novelty of the 
results in the paper lies in that 

(1) We give the first analysis for the general second order elliptic eigenvalue problem; 

(2) We remove the uniformity condition on the mesh used in the literature; 

(3) We prove the first result for the case where the eigenfunction is singular in the sense 
u G H l D {ti) n H 1+S (tt) with < s < 1; 

(4) We show the first result for the mixed boundary condition with \Tn\ ^ 0; 

(5) For both two and three dimensional cases, we only give the analysis for the rect- 
angular mesh. It is straightforward to see that the analysis is valid for the general 
quadrilateral mesh in two dimensions and the general hexahedron mesh in the three 
dimensions. 

5.2. The nonconforming rotated Qi element in any dimension. We only consider 
the case where the bilinear form a(u,v) = (Vu, Vt^arm- However, the mixed boundary 
condition is permitted. Denote by Qrq(K) the nonconforming rotated Q\ element space on 
the element K e % which reads [TjJ [79] 

(5.19) Qrq(K) := Pi(K) + span{xi - xj,xj - x\, ■ ■ ■ ,x 2 n _ x - x 2 n }, 

with the space Pe{K) of polynomials of degree < £ over the element K. For any v G H 1 ^), 
we define the following face functional 

(5.20) T f (y) = ±-jvdf 

where / is the n — 1 dimensional hyper-plane with the diameter hf. The nonconforming 
rotated Qi element space Vh is then defined by 

Vh ■= {v G L 2 (f2) : v\k G Qrq(K) for each K G Th, v is continuous with respect 

to Tf for all internal n — 1 dimensional hyper-planes / , and J~f{v) = for all / on T 

For the nonconforming rotated Q\ element, we define the interpolation operator : 
H l D {Sl) V h by 

(5.21) J Il h vdf = vdf for any v G / G H h . 
A similar argument of Lemma 15.51 below shows that 
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Lemma 5.2. There holds 

(5.22) \\u — H h u\\ L 2( K j < h 1+s \u\ H x+^ K )for any u G H l+S (K) with < s < 1 andK eT h - 

Lemma 5.3. For the nonconforming rotated Qi element, it holds the hypothesis (H4) when 
u G Hj-f Pi H 1+S (tt) with < 3 < 1. 

Proof. Since A h Vh = with the operator A h defined elementwise, we use the integration by 
parts to prove that ah{u — UhU,Vh) = for any Vh G Vh- Then, the desired result follows 
immediately from Lemma 15.21 with k = m = j = l, s = s, As = 1 — s, S = s, and 
AS = 1. □ 

In the appendix C, we show that there exist meshes such that h s < \\u — Uh\\h holds when 
u G H 1+S (Q) with < 5 < 1. Then, we have that the result in Theorem 14.21 holds for this 
class of elements. 

The theoretical result was demonstrated by the numerical results reported in [68J: the 
discrete eigenvalue will be smaller than the exact one for the singular case while the discrete 
eigenvalue will possibly bigger than the exact one for the smooth case. 

Remark 5.4. The extension of the analysis and results herein to the Crouzeix-Raviart 
element [251 179] in any dimension is straightforward. We refer the interested readers to [96J 
for further details. 

5.3. The enriched nonconforming rotated Q\ element in any dimension. Here we 
restrict ourselves to the case where the bilinear form a(u,v) = (Vit, Vi> )l 2 (u)- However 
the mixed boundary condition with |rv| ^ can be equally imposed as well. Denote by 
Qeq{K) the enriched nonconforming rotated Qi element space defined by [M] 

(5.23) Qeq(K) = Qrq(K) + span{a;i + xj, xj + xl,--- , x\_ x + x 2 n }. 
The enriched nonconforming rotated Q\ element space Vh is then defined by 

Vh ■= {v G L 2 {VL) : v\k G Qeq(K) for each K G Th,v is continuous with respect 

to J-f for all internal n — 1 dimensional hyper-planes / , and J~f(v) = for all / on T 

For the enriched nonconforming rotated Q\ element, we define the interpolation operator 
Ii h : Hb(tt) -> V h by 

/ Tlhvdf = / vdf for any v G H^Q), f G H h , 

(5.24) J ' Jf 

/ Tlhvdx = / vdx for any K ETh- 
Jk Jk 

For this interpolation operator, we have 
Lemma 5.5. There holds that 

(5.25) \\u — Hhu\\L 2 (K) ^ h 2 \u\H 2 (K)for any u G H 2 (K) and K G Th, 

(5.26) ||w — n ft n|| L 2( X) < h l+B \u\ H i+s( K )for any u G H 1+S (K) with < s < 1 and K G Th- 
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Proof. Since u — IT^m has vanishing mean on each edge of K, it follows from the Poincare 
inequality that 

\\u - U h u\\ L 2 {K) < h K \\V(u - U h u)\\ L 2 {K) . 

Then the desired result follows from the usual interpolation theory and the interpolation 
space theory for the singular case u G H 1+S (K). □ 

Lemma 5.6. For the enriched nonconforming rotated Q\ element, it holds the hypothesis 
(HA). 



Proof. We define the space Qk 



( a u + a 12 Xi \ 

«21 + a 2 2^2 



with free parameters an, a 12, • ■ ■ , a„i, a„2. 



\ a n i + a n 2X n ) 
From the definition of the operator LT^, we have 



(5.27) (V(tt - Tl h u), iI>)l*{k) = 0, for any if> G Q K . 

Let Vh be the piecewise gradient operator which is defined element by element. Since 
Vh^-hu\x £ Qk, this leads to 

(5.28) (y h u h u)\ K = p K (yu\ K ), 

with the L 2 projection operator P K from L 2 (K) onto Qk- This proves a,h(u — U^u, Uh) = 0. 
It remains to show the estimate in (H4). To this end, let n be the piecewise constant 
projection operator. Without loss of generality, we assume that p is piecewise constant. It 
follows from the definition of the interpolation operator LI^ that 

||p 1/2 n fe M||^ 2(n) - \\p 1/2 u\\ 2 L 2 {Q) = (p(U h u - u),U h u + u) L 2 (u) 

(5.29) = (p(n h tt - u), n h u + u- n (n^w + u)) L2{n) 

< h\\p 1/2 (Tl h u - u)\\ L 2 { n)\\V h (Tl h u + u)\\ L 2 {n) , 

which completes the proof of (H4) with k = j = 2, m = 1, s = S = 0, As = AS = 1, for 
the case u G H l D (Vt) n H 2 {Vl); with k = j = m = 1, s = S = s, As = 2 — 5, AS = 1, for the 
case u G H l D (Q) n H 1+S (tt) with < s < 1. □ 

In the appendixes A and C, we achieve that h < ||it — Uh\\h when u G H 2 (Q) and that 
there exist meshes such that h s < \\u — Uh\\h holds when u G H 1+S (Q) with < s < 1. 
Therefore, we have that the result in Theorem 14.21 holds for this class of elements. 

Remark 5.7. The enriched nonconforming rotated Qi element method of the eigenvalue 
problem in two dimensions with n = 2 has already been analyzed in [58J. 

6. The Adini element in any dimension for the fourth order operators 

Let the boundary dfl be divided into two parts: T D and T s with |Fd| > 0, and T D UT S = 
dfl, and T D fl T s = 0. Define the space 

V := {v G H 2 (Q),v = — = on T D , v = on T s }, 

ov 
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with v the unit normal vector to dQ. The fourth order elliptic eigenvalue problem reads: 
Find (X,u) eRxV, such that 

(6.1) a(u,v) = \(pu,v)i2(n) for any v G V and j| > o 1/ ' 2 w || z, 2 (ri) = 1 • 

We assume that the bilinear form a(-, -) satisfies the conditions of (12.21) and that p G L°°(Q) 
is a positive function. 

Let J~h be a regular n-rectangular triangulation of the domain Q C R n with 2 < n and 
denote by QAd{K) the nonconforming Adini element space [79] defined by 

(6.2) QAd(*0 = Qi{K) + ^n{x 2 i q,qeQ l {K),i = !,-■■ ,n}. 

The nonconforming Adini element space Vh is then defined by 

Vh := {f G L 2 (fi) : t>|x G Q^-fO for each K G 7/i,f and Vt> is continuous 
at the internal nodes, and vanishes at the boundary nodes on Tp, 
v vanishes at the boundary nodes on T^} . 

The finite element approximation of Problem (16.11) reads: Find (A^, Uh) G R x Vh, such that 

(6.3) a h (u h ,Vh)L2{n) = Xh(pu h ,v h ) L 2 {n) for any v h G V h and \\p 1/2 u h \\ L 2^ = 1 , 

where a h (-, •) is the discrete counterpart of a(-, •) which is defined element by element. 

For this class of elements, we have (H1)-(H3), see [79] for the details. In particular, by 
the abstract theory in Theorems 13. 2[ 13.41 and 13.51 we have 

(6.4) | A - Afc| + \\p 1/2 (u - u h )\\ L 2 {n) + h\\u - u h \\h < h 2 , 

provided that u G V R H 3 (Q). 

In order to show the hypothesis (H4), we let Ph be the Ritz-Galerkin projection operator 
defined in (13.231) . Then, we have 

(6.5) \\u-P h u\\ h < inf \\u - v h \\ h < h 1+s , 

provided that u G V fl iJ 3+3 (f2) with < s < 1. In order to derive the L 2 error estimate of 
u — PhU, we introduce the follow auxiliary problem: Find w & V such that 

(6.6) a(w,v) = {p{u - P h u),v) L 2 {n) for any v G V . 
We assume that the above dual problem admits the following regularity 

(6-7) IIHI-ff 3 (fi) ~ \\P 1/2 ( U ~ p hu)\\ L 2 {n) . 

Then we have 

(p(u - P h u),u - P h u) L 2 (n) 

(6.8) 

= (p(u - P h u),u - P h u) L 2 {n) - a h (w, u - P h u) + a h (w, u - P h u), 
which leads to 

(6.9) \\p l l 2 {u-P h u)\\ L 2 {n) <h 2+ s. 
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Lemma 6.1. For the Adini element defined in (16. 3ft . there holds the hypothesis (H4) such 
that 

(6.10) \\p 1/2 (u - P h u)\\ L 2 {n) < h 2+s ifuEVf) H 3+S (ty with < s < 1 . 

Proof. Let 11^ = P h in (H4). Since ah(u — Ph,v h ) = for any Vh E Vh, it follows from (16. 9 j) 
that (H4) holds with jfe = j = 3, m = 2, <S = 5, s = 0, < As = S, AS = 1. □ 

Then, the unifying theory in Theorem 14.21 proves the following lower bound. 

Theorem 6.2. Let (X,u) and (Xh,Uh) be the solutions of problems (16.11) and ( 16.31) . respec- 
tively. Assume that h 2 < \u — Uh\ 2 h and u E V fl H 3+S (Q) with < s < 1. JTien, 

(6.11) A ft < A, 
provided that h is small enough. 

We observe from our extensive numerical results jH] that the eigenvalue produced by the 
Adini element will be smaller than the exact one provided that its associated eigenfunction is 
smooth enough in the sense that u G H 3+S (Q) with s > 0; otherwise, the discrete eigenvalue 
will possibly approximate from both above and below to the exact one. This explains the 
numerical results given by [73] where some discrete eigenvalues are smaller than the exact 
ones and the others are greater than the exact ones. We numerically test by some high order 
methods the regularity of the cases where the discrete eigenvalues are greater than the exact 
ones and find that the corresponding eigenfunctions are singular in sense that u G H 2+S (Q) 
with s < 1. 

7. The Morley-Wang-Xu elements for the 2m-TH order operators 

This section studies the 2m-th order elliptic eigenvalue problem defined over the bounded 
domain f2 C lR n with 1 < n and m < n. Let ) be the multi-index with 

n 

\k\ = ^2 Ki, we define the space 

8=1 

(7.1) V = {v G L 2 (n), ^ G L 2 (n), \k\ < m,v\ dn = ^\ dn = 0, i = 1, • • • , m - 1}, 

with v the unit normal vector to dVl. The partial derivatives |^ are defined as 

qk v q\k\ v 

dx K dx* 1 ■ ■ ■ dx%™ 

Let D e v denote the m-th order tensor of all £-th order derivatives of v, for instance, I = 1 
the gradient, and i = 2 the Hessian matrix. Let C be a positive definite operator with the 
same symmetry as D m v, the bilinear form a(u,v) reads 

(7.3) a(u, v) = (a, D m v) L 2 (n) and a = CD m u, 
which gives rise to the energy norm 

(7.4) IMly := a ( M ) u ) f° r an y u e v , 
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which is equivalent to the usual \u\h™(q) norm for any u G V . 

The 2m-th order elliptic eigenvalue problem reads: Find (A, u) G R x V with 

(7.5) a(u,v) = \(pu,v) L 2 (n) for any v G V and \\p 1/2 u\\ L 2 {n) = 1, 

with some positive function p G L°°(Q). 

We consider the Morley-Wang-Xu element in [87] and apply them to the eigenvalue prob- 
lem under consideration. Let Th be some shape regular decomposition into n-simplex of 
the domain Q. Denote by H n -i t h, i = 1, • • • , n, all n-i dimensional subsimplexes of Th with 
v n -ij any unit normal vector to / G H n -i,h- Let [•] denote the jump over /. For any n-i 
dimensional boundary sub-simplex /, the jump [•] denotes the trace restricted to /. As 
usual, Kk is the mesh size of K G Th, and hf the diameter of / G r H. n -i,h- Given K G Th, 
let dK denote the boundary of K. The Morley-Wang-Xu element space is defined in [87J, 
which reads 

(7.6) V m , h ■= {v G L\Q),v\ K G P m (K), / [— ^-]df = 0,V/ G = V" ,™}- 

Define the discrete stress ah = CDJ^Uh, the broken versions dh{-, •) and || • \\c h follow 

ah(uh, v h ) ■ = {&h, D^v^mn), for any u h , v h G V + V^ ft , 

||tt h ||h : = a h (u h , u h ) for any u h G V + V^j, , 

where is defined elementwise with respect to the triangulation Th- 

The discrete eigenvalue problem reads: Find (Xh,Uh) GRx V m ^, such that 

(7.7) a h (u h ,v h ) = \h(pu h ,v h )L2(n) for any v h G V m>h and ||p 1/2 W/,||ra(n) = 1- 

The associated interpolation operator for the space V m ,h is defined by: Given any v G V, 
the interpolation il m t> G is defined as 

(7-8) / m _" rf/ = / - d/, for any / G W n -»,fc , i = 1 , • • ■ , m . 

Jf ov n _ i f J f av n _ i f 

For this interpolation, we have the following approximation 

(7.9) \\p 1/2 (u - U m u)\\ L 2 (n) < h m+B \u\ Hm +s (n ) for any u G V n H m+s (Vl) with < s < 1 . 

It is straightforward to see that the hypothesis (H1)-(H3) hold for this case [79j[87]. Then, it 
follows from the general theory for the error estimate of the eigenvalue problem in Theorems 
13.21 13.41 and 13.51 and the properties of the nonconforming methods used herein [79| 187] that 

(7.10) \\u - u h \\ h < h 5 and ||u - u h \\ L ^(n) < h 2s , 
provided that the eigenfunction u G V D H m+S (Q) with < s < 1. 

Theorem 7.1. Let (X,u) and (Xh,Uh) be the solutions of problems (17.51) and (17. 7p . respec- 
tively. Then, 

(7.11) X h < A, 
provided that h is small enough. 
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Proof. It follows from the definition of the interpolation U m in (17 .8p that ah{u — U m u, Vh) = 
for any Vh G Vh- The hypothesis (H4) follows immediately from ( 17. 9ft . In addition, in the 
appendixes A and C, we show that h < \u — Uh\h when u G V D if m+1 (f2) and that there 
exist meshes such that /i s < ||m — w^lU holds when w e i? m+s (f2) with < s < 1. Then, the 
desired result follows from Theorem hL2J for m > 2. □ 



8. New nonconforming elements 

In this section, we shall follow the condition (H4) and the saturation condition in Theorem 
14.21 to propose new nonconforming finite elements for the second order operator. This is 
of two folds, one is to modify some nonconforming elements in the literature such that the 
modified one will meet the condition (H4), the other is to construct new nonconforming 
elements. 

8.1. The enriched Crouzeix-Raviart element. To fix the idea, we only consider the 
case where n = 2 and note that the results can be generalized to any dimension. Let Th 
be some shape regular decomposition into triangles of the polygonal domain Q C 1R 2 . Here 
we restrict ourselves to the case where the bilinear form a(u,v) = (Vw, Vu)i, 2 (fi) with the 
mixed boundary condition |Tjv| ^ 0. 

Note that the original Crouzeix-Raviart element can only guarantee theoretically the lower 
bound of the eigenvalue for the singular case in the sense that u G H 1+S (Q) with < s < 1. 
To produce the lower bound of the eigenvalue for both the singular case u G H 1+S (Q) and 
the smooth case u G H 2 (Q), we propose to enrich the shape function space by x\ + x\ on 
each element. This leads to the following shape function space 

(8.1) Qecr(K) = Pi{K) + span{x 2 + x\} for any K G T h . 

With the functional J 7 / from ( 15.201) . the enriched Crouzeix-Raviart element space Vh is then 
defined by 

Vh := {y G L 2 (Q) : v\k G Qecr(K) for each K G Th,v is continuous with respect 
to J "f for all internal edges / , and J~f(v) = for all edges / on Y D j . 

For the enriched Crouzeix-Raviart element, we define the interpolation operator 11^ : Hq(Q) — > 
V h by 



5.2) 



J Hhvdf = J vdf for any v G if^(O) for any edge / 



Uhvdx = / vdx for any K G Th- 
k J K 



For this interpolation operator, a similar argument of Lemma 15.51 leads to: 
Lemma 8.1. There holds that 

(8.3) ||w - Il h u\\ L 2 {K) < h 2 \u\ H 2 [K) for any u G H 2 (K) and K E T h , 

(8.4) \\u - U h u\\ L 2( K ) < h 1+s \u\ H i+s( K )for any u e H 1+B (K) with < s < 1 and K e T h - 
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Lemma 8.2. For the enriched Crouzeix-Raviart element, it holds the hypothesis (H4)- 

Proof. We follow the idea in Lemma 15.61 to define the space Q K = 11 ' 1Z 1 with 

\ a 2 i + a 12 x 2 J 

free parameters an, a 2 i, a\ 2 - From the definition of the operator IL^, we have 

(8.5) (V(it - U h u), ip) L *(K) = 0, for any G Q K . 
Indeed, we integrate by parts to get 

(V(tt - U h u), ^)ls(k) = ~{u - U h u, div tp) L 2 {K) + ^ ( u - ^-hu)il> ■ v f ds. 

fCdK J f 

Since divi/> and if) ■ Vf (on each edge /) are constant, then (18.51) follows from (18.21) . Since 
VhJlhu\K G Qk, the identity (18. 5p leads to 

(8.6) {V h U h u)\ K = P K {Vu\ K ), 

with the L 2 projection operator Pk from L 2 (K) onto Qk- This proves ah{u— ILjU, u^) = 0. It 
remains to show the estimate in (H4). Then, it follows from the definition of the interpolation 
operator that 

, g ^ \\p 1/2ll hu\\ 2 L 2 in) - \\p 1/2 u\\ 2 L 2 {Q) = (p(Jl h u - u), U h u + u) L 2 (n) 

= (p(IL h u-u), U h u + u- IIo(n fc ti + u)) L 2( Q) , 

which completes the proof of (H4) with k = j = 2, m = 1, s = S = 0, As = AS = 1, for 
the case u G Hp(Q) fl H 2 (tt); with k = j = m = 1, s = S = s, As = 2 — 5, AS = 1, for the 
case u G H^Q) n H 1+S (Q) with < s < 1. □ 

In the appendixes A and C, we have proven that h < \\u — Uh\\h when u G H 2 (Q) and 
that there exist meshes such that h s < \\u — Uh\\h holds when u G H 1+S (Q) with < s < 1. 
Whence, we have that the result in Theorem 14.21 holds for this class of elements. 



8.2. A new first order nonconforming element. With the condition from Theorem 14. 21 
a systematic method obtaining the lower bounds for eigenvalues is to design nonconform- 
ing finite element spaces with good local approximation property but not so good global 
continuity property. To make the idea clearer, we propose a new nonconforming element 
that admits lower bounds for eigenvalues. Let l~h be some shape regular decomposition into 
triangles of the polygonal domain fid 2 . We define 

Vh := [v G L 2 (Q) : v\k £ ^(-^O for each K G 7h,v is continuous with respect 

to Tf for all internal edges / , and Tf{v) = for all edges / on T D } . 

For this nonconforming finite element space, one can follow the same line for the Wilson 
element to show that the discrete eigenvalue by it for a general elliptic operator is smaller 
than the exact one provided that the eigenfunction u G V fl H 1+S (Q) with < s < 1 or 
u G V n H 2+S (tt) with s > 0. 
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Next, we give an example to show that the lower bound is not promised theoretically. 
Define 

V h := {v G L 2 (Vl) : v\ K G P X {K) © span{3(L? + L\ + L|) - 2} for each K G 7h, 
t> continuous with respect to J 7 / for all internal edges / , and J~f{v) = 
for all edges / on Td} , 

where L\, L 2 and L 3 are barycentric coordinates of the element K. This nonconforming 
space is proposed by Brezzi and Marini [20]. Although it looks very like the space defined 
in (18.21) . it is generally impossible to prove that it admits the lower bound for the eigenvalue 
when the eigenfunction is smooth. The reason lies in that the continuity of the enriched 
nonconforming local bubble function 3(Lf + L\ + L\) — 2 is of order 2 over each internal 
edge and therefore better than the local approximation property. 

These examples conclude that we can derive the lower bound for the eigenvalue through 
the local enrichment. But not all enrichments work and some well-tailored balance between 
the global continuity and the local approximation is obviously dispensable. Such a balance 
is nothing else but the condition that the local approximation property is better than the 
global continuity property. 

We end this section by giving some remarks. We note that both the original Crouzeix- 
Raviart element and the nonconforming rotated Qi element can not assure theoretically 
the lower bounds for eigenvalues when the corresponding eigenfunctions are smooth. Our 
comment is that their local degrees of freedom are not enough when compared to the global 
continuity. Hence we enrich them by some local nonconforming bubble functions on each 
element and find that their enriched versions produce the lower bounds for both the smooth 
case and the singular case. The analysis in a recent paper [97] indicates that such an 
enrichment is in some sense necessary for the Crouzeix-Raviart element. In fact, the result 
therein says that A^ < A provided that ||w — Uh\\h > CXh with C > l/\^12 being a positive 
constant. On the contrary, the Crouzeix-Raviart element will possibly give the upper bound 
for the eigenvalue. The similar argument equally applies to the nonconforming rotated 
Qi element. These comments are actually confirmed by the numerical examples for the 
nonconforming rotated Q\ element [681 EH]- 

9. Conclusions 

In this paper, we propose a systematic method that can produce the lower bounds for 
eigenvalues of elliptic operators. With this method, to obtain the lower bound is to design 
some nonconforming finite element space with enough local degrees of freedom when com- 
pared to the global continuity. We check that the most used nonconforming methods in 
the literature possess this promising property. We also propose some new nonconforming 
methods with this feature. In addition, we study systematically the saturation condition 
for both the conforming and nonconforming methods. 
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Appendix A. The saturation condition where the approximation error is 

DOMINANT 

In the following three sections, we shall prove, for some cases, the saturation condition 
which is used in Theorem 14.21 The error basically consists of two parts: the approximation 
error and the consistency error due to the nonconformity of the discrete space. Therefore, in 
this section, we analyze the case where the approximation error is dominant; in the following 
section, we are concerned with the case where the consistency error is dominant; in appendix 
C, we give some comments for the case where the eigenfunction is singular. 



A.l. An abstract theory. Let u G V D H m (Q) be the eigenfunction of some 2m-th order 
elliptic operator. Let Vh be some fc-th order conforming or nonconforming approximation 
space to H m (Q) over the mesh 7~h in the following sense: 

(A.l) sup inf ||-D^(w — Vh)\\L 2 (n) ^ h k \v\ H m+k for some positive integer k. 

V (zH m + k (n)nv VheVh 

Then the following four conditions are sufficient for the saturation condition: 

Jl. There exists a local interpolation operator H m +k from the space H m+k (Q) to some 
high order finite element space than the finite element space Vh under consideration; 
J2. The following Poincare inequality 

(A.2) \\D%(u - U m+k u)\\ L 2 {n) < h k \\D™ +k (u - Il m+k u)\\ LHnh 

holds for the local interpolation operator U m+ h] 
J3. The following basic approximation property 

(A.3) \\Dt + \u - U m+k u) || L2(n) when h->0, 

holds for the local interpolation operator U m+ h] 
J4. At least one fixed component of D^ +k Vh vanishes for all Vh G Vh and while the L? 
norm of the same component of D m+k u is nonzero. 

We recall that D e v denote the m-th order tensor of all £-th order derivatives of v, for 
instance, £ = 1 the gradient, and £ = 2 the Hessian matrix, and that D l h are the piecewise 
counterparts of D e defined element by element. 

Theorem A.l. Under conditions J1-J4, there holds the following saturation condition: 
(A.4) h k <\\D%(u-u h )\\ Lm . 

Proof. By the condition J4, we let 9T denote the multi- index set such that \k\ = m + k for 
any k G 9t and that 

(A.5) = for any K G % and v h G V h while \\-^\\ L 2 {n) ^ 0. 
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Whence it follows from the triangle inequality and the piecewise inverse estimate that 

E u d K u 2 \- \- u d K (u - u h ) ||2 

"Q x K"L?(n) ~ 2^ 1^ II q x k Wl*{k) 

(A.6) , (^{u - n m+fc M) 2 9 K (n m+fc M - 2 x 

- 1^ I" <^ Hi 2 W + I' Q^K \\L\K)) 

KemKeT h 

< \\D™+ k (u - U m+k u)\\ 2 L2{n) + h- 2k \\D-(U m+k u - u h )\\ 2 L2(U y 
By the Poincare inequality in the condition J2 and the triangle inequality, we obtain 

(A.7) E £ W D T + % - n m+k u)\\ 2 LHn) + h-' k \\D™(u - u h )\\l 2(n) . 

new 

Finally it follows from the condition J3 that 

QK, 

^ dx 
K£<n 

when the meshsize is small enough, which completes the proof. □ 



(A.8) 



A. 2. The Morley-Wang-Xu element. We first consider the case where m < n — 1. 
Assume that the exact eigenfunction u G H m+1 (Q). We let U m+ iu denote the canonical 
interpolation defined as in (17.81) in the discrete space V m+ i t h from (17.61) . Integrating by parts, 
we get from the definition of H m +i that 

(A.9) D^ +1 U m+1 u = U D m+1 u, 

we recall that Ho is the piecewise constant projection operator with respect to the triangu- 
lation Th of the domain Q C M. n , see [501 Lemma 3.9] for more details. Since fjD™(u — 
U m+ iu)df = for any n — 1 dimensional sub-simplex /, we have the following Poincare 
inequality. 

(A.10) \\Dl\u - U m+1 u)\\ L 2 (u) < h\\D™ + \u - IWnOHi^n). 

This establishes the conditions Jl and J2 for this class of elements with the integer k = 1. 
Now it follows from (1A.9|) that 

(A.ll) \\D^ +1 (u-U m+l u)\\ L 2 {u) = \\D m+1 u-Tl D m+1 u\\ L2m . 

From the density of piecewise constant functions in L 2 (Q) we have that \\D m+1 u— n -D m+1 -u||i2(^) 
tends to zero as h goes to zero, which proves the condition J3. 
For this class of elements, we have 



(A.12) \\D™ +1 v h \\ L 2 {n) = for any v h G V m , h . 

We observe that ||-D m+1 w||£2(Q) can not be equal to zero since u is the eigenfunction of the 
differential operator A. Otherwise u will be some polynomial, which is impossible to be the 
eigenfunction of any differential operator. This ends the proof of the condition J4. 

Remark A. 2. The above analysis applies equally to the conforming linear element for the 
second order problem. 
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Next we consider the case m = n. For this case, there does not exist the m + 1st Morley- 
Wang-Xu element in the literature [87]. Therefore we need to construct a new interpolation 
such that the commutative property ( 1A.9I) holds. For ease of presentation, we only consider 



the case where n = 2. Given any element K with three vertexes V\, Vi y and V3. Let e±, e2, 
and e3 be three edges, which are numbered counterclockwise and are associated to the unit 
normal vectors i — 1 , 2 , 3, respectively. For any v G H 3 (K), we define Hrv G P%{K) by 

(A.13) ggg^ | P< = 0,i = l,2,€=l,2,3, and f V) ds = 0, 1 = 1 , 2 , 3, 

and (n^t>)(xo, yo) = ^(^0,2/0) where (xo,2/o) is the barycenter of K. It is straightforward to 
see that the interpolation is well-defined and that the commutative property ( 1A.9I) holds. 



Then we can use the same argument for the case m < n — 1 to prove the conditions J1-J4 
with k = 1 and Il m+ i|^: = IIk for any K G l~h- 



A. 3. The enriched Crouzeix-Raviart element. For this element, we have not the prop- 
erty || V^u^H^m) = where V 2 h Uh is the piecewise Hessian matrix of Uh defined by elemen- 
twise. For ease of presentation, we only consider the case where Q = [0, l] 2 and u G (Q). 
Suppose u € H 2 (Q) and let n 2 M be the interpolation of u by the Morley element in the 
two dimensions. A similar argument for the Morley- Wang-Xu element shows the conditions 
J1-J3 for this element with k — 1. 

Let di2,h denote the piecewise counterpart of the differential operator ^f^. We have 
d\2,hVh = for any Vh € Vh- To finish the proof of the condition J4, we need to show that 
||^^IU 2 (f7) can n °t vanish. Otherwise, u should be of the form u(x, y) = f(x) + g(y), where 
f(x) is some function of the variable x and g(y) is some function of the variable y. Now 
the homogenous boundary condition concludes that f(x) = C\ and g(y) = C2 for some 
constants C\ and C2, which are impossible to be eigenfunctions of any differential operator. 
Therefore the condition J4 holds for this case. 

The same argument applies to the nonconforming Qi element, the enriched nonconforming 
rotated Q\ element, and the conforming Q\ element in any dimension. For these cases, we 
have to use the interpolation operator associated to the n-rectangle Morley element, whose 
degrees of freedom are the value of the normal derivative at the centric point of each n — 1 
dimensional face and the function value at each vertex [SB]. For each element K, the 
shape function space is Qi(K) + spanjxf , x 2 ,,--- , x 2 w x±, x\, ■ • • , x\}. For this case, we have 
not the commutative property (1A.9jl . However, we can use the approximation property of 
the n-rectangle Morley element space and the density argument to finish the proof of the 
conditions J1-J4. 

A. 4. A special proof for the case m = 1. For this case, we shall use the bubble function 
technique to prove the saturation condition. The bubble function technique is frequently 
employed in the community of a posteriori error analysis [OH [T]. For a posteriori error 
analysis of the eigenvalue problem, we refer to the references [321 ES, HHJ ETJ [59] . 

Given any K G we let bx G Hq(K) denote the usual bubble function, which is a 
polynomial and equal to 1 at the barycenter of K. We define wk = bxXhUh, consequently 
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have wk G Hq(K) and 

(A.14) \w\hi(Q) < h^WXhllhWtf^y 

Whence, we derive, by using the elementwise inverse estimate, the condition A^-u^ = (we 
recall that is the piecewise counterpart of the Laplace operator A, which is defined 
element by element) and integrating by part, 

W^hUh\\ 2 L 2 {K) {\hU h ,W K )L?(n) 

= (X h u h - \u + \u - Am + Am + A h u h , w K ) L 2 {K) 

(A.15) 

< \\\ h u h - \u\\ L 2 {K) \\w K \\ L 2 {K) + a h (u - u h ,w K ) 

< \\\ h u h - \u\\ L 2 {K) \\w K \\ L 2 {K) + h^\\u - u h \\ K \\w K \\ L 2 {K) , 

where \\u — Uh\\ic is the restriction of \\u — Uh\\h over K. Note that 1 1 A^Uft, — \u\\l2^ K ) is of 
high order provided that u G H 1+S (£l) with < s. A summation of all the elements yields 

(A.16) h\\X h u h \\ L 2 {n) < \\u - u h \\ h . 

Note that we do not need the condition u G H 2 (Q) for this analysis. 

Appendix B. The saturation condition where the consistency error is 

DOMINANT 

In this section, we prove the saturation condition for the case where the consistency error 
is dominant. As usual it is very implicated to give an abstract estimate for the consistency 
error in a unifying way. Therefore, for ease of presentation, we shall only consider separately 
the nonconforming Wilson element and the new first order nonconforming element proposed 
in this paper. However, the idea can be extended to other nonconforming finite element 
method. 

B.l. The Wilson element. For each K G we introduce the following affine invertible 
transformation 

F K : K ->• K,Xi = h x . tK £i + 
with the center (x^x®) and the lengths 2h Xii K of K in the directions of the rcj-axis, and 
the reference element K = [—1, l] 2 . Take v h \ K = a ltK (l — £ 2 ) + a 2 ,x(l — £i) w ith some 



parameters a^K, we have 

(B.l) HVtfcHi^ = ^Kk^ + alJ^). 

For any wk G P2(K), we choose a^K = 9 dx f and a,2,K = 9 dx i c ■ A direct computation gives 

(B.2) J -7^v h ds = -(a\ K + a\ K )h xl:K h X2tK . 

Therefore 

/T-, n\ f ® U f d( U ~~ W k) 8/2 2 r n / 

(B.3) / — v h ds = / v h ds + -{a 1K + a 2 K )h xuK h X2iK for any w K G P 2 {K ). 

JdK av JdK av <J 
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A summation over all elements yields 



E IdK^ Z h ds E (§( fl l,y + a 2.K) h xuKh X2 ,K + j dK {U d ™ K > V h ds) 

(B.4) sup ||V7 - |j > — 575- , 

y ( a l K^2,K /h Xl ,K + a 2 , K h xlj K/h X2j K)) 

KeT h 

for any % G iV-^O and = ai,jc(l-^i) + a2,A'(l-^i)- For u G iT 2+s (n) with < s < 1, 



e Pa' 
this proves that 

rr> r\ KTi . „ , ,,,£7 M 2 ..OIL.. vl/2 „ , i +sl 1 

(B.5) sup ^c^fli || +|| || L2(n) ) -C 2 fc + M^+ S( n), 

for some positive constants C\ and C2. 

Since the error || Vh(u — Uh)\\L 2 (n) consists of both the approximation error and the consis- 
tency error, the lower bound just established for the consistency error proves the saturation 
condition when the meshsize h is small enough provided that u G H 2+S (Q) with < s < 1. 
In fact, the analysis proves that 

(B.6) h< \\u-P' h u\\ h , 

where the projection operator P' h is defined as in (13 .4ft . Note that 

(B.7) \\P' h u - u h \\ h < \\Xu - \ h u h \\ L 2 {n) < h 2 , 

provided that u G H 2 (Q). This gives 

(B.8) /i< ||u-u h || A . 

The analysis herein weakens the regularity requirement on u from [23], where u G W 3,OQ (Q) 
is needed. 

B.2. The new nonconforming element. Given any element K, let ei, e 2 , and ej, be three 
edges, which are numbered counterclockwise and are associated to the unit normal vectors 
Ui, i — 1 , 2 , 3, respectively. We denote by (fi,K, i = 1,2,3, the interior basis functions of 
•PbC^O such that 



(B.9) 



/ (fj )K ds = 0,i,j = 1,2,3. 

J e, 



we define (fix = Yl Ci t K<fi,K, for some interpolation coefficients c^k- Let Ck = {ci,k, C2,K, c 3,x) T , 
i=i 

we shall prove that there exists some interpolation n 2 M G P^{K) such that 

(B.10) \\V{U 2 u-u)\\v {K) < h l+s \u\ m+ s {K) 

provided that u G H 2+S (K) with < s < 1, and that we have the following estimate 



3 ,2 



(B.H) sup 3 \ ' ' / > C7/i||V 2 n 2 n|| L2W , 

O^CkSK 3 II V0A-||z2(K) 
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where Tj = (—Vj,2, is the unit tangential vector to the edge e r To this end, we define 
the matrix Ak G M 3x3 by the components 

(B.12) dij = / (p j>K sds,i,j = 1,2,3. 



Note that a^j = 0(h 2 ). We recall that Li, L 2 and L 3 are barycentric coordinates of the 
element K. We know that bx = 3(Lf + L\ + L|) — 2 is the unique bubble function (up to 
some multiplier constant) such that 

(B.13) \ b K s l ds = 0,j = 1,2, 3, £ = 0,1. 

Whence the rank of the matrix A K is 2. Denote £ x = (B 1>K , B 2>K , B 3 , K ) T = (g^, g^, gf) T 
Then the estimate (IB. Ill) can be rewritten as 

(B.14) sup ,l\ A f K > Ch\\V 2 U 2 u\\L H K). 

To achieve both flRTUl) and (lB~14j) for some Ii 2 u G P 2 (iT), we first let n 2 u G P 2 (iT) be 
some interpolation with (1B.10I) . If A K B K ^ ( B K is defined in the same way as B K ), 



then there exists some Ck G M 3 such that *4.xCx = Bk and therefore we take U 2 u = n 2 w 
and i?^ = B K ; otherwise, we slightly modify Yl 2 u to get n 2 w such that A K B K ^ and 
||V(n 2 w — n 2 ti) || z, 2 (^) ^ h 1+s \u\H 2 +s(K)- in both cases, we can choose the vector Ck such 
that 

3 3 

(B.15) !|V0dl V) < E l c ^| 2 S ^ E l^l 2 - 



i=l i=l 



Therefore 



(B.i6) sup ^ .^*f K > ch 2 (J2 l^l 2 ) 1/2 « C%|| v 2 n 2M || i2W . 



(VCVeM 3 ||V0x||L2(iC) 



8=1 



Finally, we have 
(B.17) 

e j^i^ e j^i^ 

SUp — — ii > > dhW V M|| L 2(n) - C 2 /l + |u|ffa4-(0), 

o^ey*, || VhWh||L2( Q ) ||Vh(. ^ifJIU 2 (n) 

provided that u G H 2+S (Q) with < s < 1 for some positive constants Ci and C 2 . Since 
II V 2 M||i2(Q) can not vanish, this proves the saturation condition. 

Appendix C. The comment for the saturation condition of the singular 

case 

We need the concept of the interpolation space. Let X, Y be a pair of normed linear 
spaces. We shall assume that Y is continuously embedded in X with Y C X and || • \\x < 
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|| • || y. For any t > 0, we define the K— functional 

(CI) K{f, t) = K{f, t, X, Y) — inf ||/ — g\\ x + t\g\ Y , 

where || • \\x is the norm on X and | • |y is a semi-norm on Y. The function K(f, .) is defined 
on M_)_ and is monotone and concave (being the pointwise infimum of linear functions). If 
< 9 < 1 and 1 < q < oo, the interpolation space (X, Y) d>q is defined as the set of all 
functions / el such that [HI ED |30] 



(C2) \f\(x, Y)e , t 



(j2 [2( s +^ ke K(f, 2-*( s + e ))]«)Vg ) o < q < oo, 

k=0 

sup 2 k ( s +^ 9 K(f,2- k ^),q = oo, 

fc>0 



is finite for some < s + e < 1. 

C.I. An abstract theory. We assume that u G H m+S (Q) with < s < 1 and Vh is some 
nonconforming or conforming approximation space to the space H m (Q). 

Then the following conditions imply in some sense the saturation condition for the singular 

case: 

J5. There exists a piecewise polynomial space V^ +lh C H m+1 (Q) such that V^ +lh C 

Kn+i h/2 w hen Th/2 is some nested conforming refinement of Th\ 
J6. There holds the following Berstein inequality 

(C.3) \v\ H ™+°+t(n) < ^ (s+e) Mtf m (n) for any v G V^ +1<h ; 

J7. There exists some quasi-interpolation operator II : V/j — > V^ +lh such that 

(C.4) \\D m (u-Uu h )\\ L2in) <h s ^ 

provided that \\D™(u — Uh)\\i^(a) < h s+e with e > and s + e < 1. 

Theorem C.l. Suppose the eigenfunction u G H m+S (Q) with < s < 1. Under conditions 
J5-J7, there exist meshes such that the following saturation condition holds 

(C5) \\D™(u-u h )\\ L2(n) . 

Proof. We assume that the saturation condition h s < ||D^(u — «fe)||L 2 (n) does not hold for 
any mesh Th with the meshsize h. In other word, we have 

(C6) \\D%(u-u h )\\ L 2 {Q) <h s+ % 

for some e > 0. In the following, we assume that s + e < 1. By the condition J7, we have 
(C.7) " % \\D'h{u - Tlu h )\\ LHn) < h s+ \ 

Take X = H m (9) and Y = H m+s+£ (n). The inequality (ICl7l) is the usual Jackson inequality 
and the inequality ( 1C.3I) is the Berstein inequality in the context of the approximation 
theory [3lJ[30]. We can follow the idea of [2U Theorem 5.1 , Chapter 7] to estimate terms 
like K(u, 2~ £<s+e )) for any positive integer £. In fact, let ^p k G V^ +1 2 - k ( s + e ) be the best 
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approximation to u in the sense that \\D m {u — ipt) ||l 2 (q) = inf || -D m (-u — t>)||L 2 (n), 

Ve m+l,2- fc ( s + £ ) 

k > 1. Let ipk = ¥k — ifik-i, k — 1, 2, • • • , where = y^o = 0. Then we have 

(C.8) p m ^|U 2 (n) < \\D m (u - <p k )\\ L i (U ) + \\D m (u - <p k -i)\\ L 2 {u) < 2- fc(s+e) . 

i 

Since <fi — ^2 ipk an d |V'o|iT m + s + e (o) = 0, it follows from (lC.7j) and (IC.3I) that 



fc=0 

(it, 2-( s+e ) £ ) < ||u - + 2- ( - s+ ^\i Pe \ Hm+ s + , 

t 

(C.9) <2-^ + 2^)^2 fc ^)||^||^ (Q) 

k=l 



(CIO) |tV Wm+ * +e( n)) e , 2 = (J][2 fc ^)^(n, 2 - fc ^))] 2 )V2 < (^[A;2 fc ^)( e - 1 )] 2 ) 1 /2. 

fc=0 fc=0 

Let 6 1 = 1 — e with e > such that e — (s + e)e > 0. This leads to 



oo 



(C.ll) | U | (ifm(n))Hm+s+ e (n))9i2 < (^[k2- k ^*°] 2 ) 1/2 < oo. 

fc=0 

This proves that u G if m+ ( 1_eo )( s+<E )(f2) which is a proper subspace of H m+s {Vl) since e — (s + 
e)eo > 0, which contradicts with the fact that we only have the regularity u 6 H m+S (Q). □ 

Remark C.2. T/ie inverse theorem was proved for the approximation by the piecewise poly- 
nomials in [92] , see also |5J. Such an inverse theorem in fact implies the above theorem for 
some special cases if we work within the Besov space instead of the Sobolev space. However, 
those results are inapplicable for the more general cases under consideration. 

C.2. The first order finite element method for the Laplace operator. We shall 
consider the first order nonconforming finite element method for the second order elliptic 
eigenvalue problem. For ease of presentation, we consider the two dimensional case with 
V = Hq(Q) and we stress the idea can be easily extended to any dimension. To this end, 
we let Vh be the conforming or nonconforming space under consideration. We assume the 
following weak continuity for any Vh G Vh'. 

(C.12) Vh is continuous on at least one point for any interior edge E or / [vh]ds = 0, 

J E 

where [vh] is the jump of Vh across the edge E. In addition we assume that 
(C.13) \\Dlvh\\ L 2 {n) = for any v h E V h . 

Let V£ be the conforming element space of H 2 (Q) over Th such that we have the inclusion: 
The shape function space on element K for Vh is some subspace of the shape function 
space on element K for V£ for any element K G Th- The existence of such a space can be 
guaranteed by Lemma IC.4I below. This establishes the conditions J5 and J6. 
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To prove the condition J7, we need to construct some quasi-interpolation operator IT : 
Vh —> V 2 . Let Vk be the restriction of V 2 on K and cpi, i = 1, , 2 • • • , £, denote its nodal 
basis functions associated with dual basis Fi, i = 1, 2 • • • , £, in the sense that F i(<fj) = 

Given v, we assume that F%(v) = 9 & v ^f l for some multi-index k 1 = (k\,k 2 ) and for some 
node p of element K. Given any node p, we define the nodal patch by 

(C.14) u p = {K eT h , P eK}. 

Given any Vh € Vh, we define the quasi-interpolation Hvh € V 2 : 
(C.15) Fi(Jlv h ) = ^Y1 f *( v >*\k), 

where p is the associated node for F % and £ p is the number of elements in uj p . We introduce 
some notation. Let Eh be the set of all edges in Th, be the set of interior edges, £d be 
the set of edges on r fl , and Sn be the set of edges on T N . 

Lemma C.3. Let the quasi-interpolation operator IT be defined as in ( 1C.15I) . Then, 
(C.16) ||V h (IIVfc — Vh)llia( ) < Yl h E\\NhVh}\\ 2 L 2 { E) f° r an V v h e V h . 

Proof. It follows from the definition of II that 

(C.17) F i(Jlv h \ K - v h \ K ) = F i(Uv h \ K ) - F i{v h \ K ) = j- ^ (F i(v h \ K ') - F i(v h \ K ))- 

The elements in u p can be rearranged as a sequence Ki, K 2 , • • , K$ with K\ = K^ p+ ± such 
that Kj and Kj+i share a common edge Ej. This yields 

(C.18) \F i{Uv h \ K ~ Vh\k)\ < 1^ i( V h\Kj - Vh\K 3+1 )\- 

3=1 

Since Fi(<fi) = 9 — 1 an d Fj{<Pi) = for Fj other than Fj, we use the scaling 

argument and the inverse estimate to obtain that 

llV^illia^ < h 2 ^. 
Since F i(Yivh\K — VhIk) = if 2 < this and (10. 181) yield 

\\Vh(Hvh-Vh)\\ 2 L 2 iK) < \Fi(u.v h \ K - v h \ K )\ 2 h 2 ^ 

|/c»|<2 

(C.19) t-l i 



M<2 j=l 



dx K * dx K% 



Next, we estimate the terms I- — h ^ K ^ _ ^ — h K ] + 1 ^ p i ^y e g^gj] use ^ e wea k continuity of 

cte K dx K 

Vh inherited from the discrete space. We first consider the case where k % = 0. Let Ej be the 
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common edge shared by Kj and Kj + i. The weak continuity of Vh across Ej shows that 

(C20) \i>h\K s (p) ~ v h \ Kj+1 (p)\ < ^ 2 ||[V h Vfc] ■ T Ej \\v*( Ej ), 

where r Ej is the unit tangential vector for edge Ej. Whence we apply the inverse estimate 
to get 

rr91 N , d Kt Vh\ Kj (p) d Kl v h \ Kj+1 (p) lKi . 1/2 

We insert the estimate (1C.21j) into ( 1C 19)) and use the fact that F i(Hvh\K — vh\k) = if 
£ p = 1 to prove 

(C.22) \\Vh(Ilv h -v h )\\l 2(K) <J2 E h E\\[VhV h ]\\h (E) , 

V E<=£ p n£n 

where S p denotes the set of edges which have p as one end point when p is a vertex of Th 
while S p denotes the edge on which £ p lies when p is not a vertex of Th- The summation 

^2 is done over all nodes on K. Then the summation over all the elements completes the 

p 

proof. □ 

Let Vh = Uh in the above lemma, we obtain 
(C23) l|V h (nu h -u A )lli»(n) £ E M[V fe ^]||£ w 

Now we shall prove that 

(C.24) h E\\[VhU h }\\ h{E) ~ - Wfe)|||a(n) + high order term. 

Such an inequality is usually referred to as the efficiency of the a posteriori error estimator 
for the nonconforming part of the error in the context of the adaptive finite element method. 
The proof can be actually found in [59] . For readers' convenience, we give a brief proof. 

We first bound the tangential part of the gradient on each edge. Let u E denote the union 
of element (s) that share the edge E and b E denote the unit edge bubble function of edge 
E in the sense that b E G Hq(uj e ), b E is equal to 1 at the mid-point of E, and that b E is a 
piecewise polynomial of degree 2. Define w E = b E \S7hU%] ■ t e , we have 

(C.25) HtVfcitfc] •TeH^) « {[V h u h ] ■t e ,w e )l^(e) = (V h (u h - u),CutIw e ) L 2 (uje) , 

where the operator Curl = (— J^, j^) T ■ The inverse estimate || Curl we\\l 2 (u e ) ~ ^e^II [Vhit/J- 
t e\\l 2 (e) shows that 

1 /2 

(C.26) h£ 1 1 [V h u h ] • r E || L 2^ E) < \\ V h (u - u h ) \\ l2( Ub ). 

To bound the normal component of the gradient on the internal edge, we set we = 
b E [VhUh] ■ v E . Integrating by parts we obtain 

^ c 2 ^ II NhU h • v E ) \\ 2 L 2 [E) « ([Vfcltfc] ■ v E , w e ) L 2 (u1e) 

= (V h (u h - u),Vw E ) L 2 (ulE ) - (div h V h {u h - u),w E ) L 2 (oJE ). 
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— 1 II 

The inverse estimate || Vwe\\l 2 (lo e ) ^ h E \\[VhUh] • ve\\l 2 (e) proves 

— 1/2 

(C.28) \{Vh{uh-u),VwE)iP(u B )\ <h E ' \\V h (u - u h )\\ L 2 { u E )\\[V h u h } ■ v e \\l2(e)- 

The second term on the right-hand side of flC.27j) can be estimated by 

| - (div h V h {u h - u), w e )l2(u b )\ 
(C.29) _ 1/2 

< h E \\\ShUh\ ■ ^e\W{e){\W h{u - u h )\\ L 2 {uiE) + || Am - X h u h \\ L 2 {uJE) ). 

Here we use the fact that (V«, V^b)^^^ = A(w, We)l 2 (w e )- We combine ( 1C.27l) -l lC.29l) to 
obtain 

1/2 

(C.30) h£ llfVftttfc] • I/js||ia(B) < ||Vfc(lt - MfcjIll^CwB) + W^ U _ ^hV>h\\V(u B )- 

From the estimates (1C.26j) and (1C.30[) . a summation of all edges proves the second inequality 

of ra . 

Finally we combine estimates (1C.23I) . (1C.24j) . and the triangle inequality to show the 
condition J7. 

C.3. The Morley-Wang-Xu element. This subsection checks the conditions J5-J7 for 
the Morley-Wang-Xu element. We need the following result from [27]. 

Lemma C.4. There exist piecewise polynomial spaces V° mh with nodal basis overTh such 
that V nt m t h are nested and conforming in the sense that V^ mh C V° mh i 2 C H m (Q) for any 
1 < n and m < n. 

This lemma actually proves the conditions J5 and J6. Let us prove the condition J7 for 
this family of elements. Let Huh G KTm+i h denote the quasi-interpolation defined as in [501 
(3.6)], cf. AC. 151) for the two dimensional case. Then we have 

(C.31) \\DZ(JLu h -u h )\\v m < \\D^(u-u h )\\ L 2 {n) + h m \\X h u h \\ L 2 (n) < h s+ \ 

We refer the interested readers to [50, Theorem 3.2] for all the details. In the proof we shall 
use the fact that ||-D^ +1 t , / l |U 2 (n) = for any Vh £ V m ^- This shows the condition J7. 

Appendix D. An economical method for the upper bound of the eigenvalue 

The idea of Section O motivates an economical method for the upper bound of the eigen- 
value. In fact, this work was finished in [3H]- We recall some details as follows. Let Uh be 
the discrete eigenvalue by the Morley-Wang-Xu element or the first order nonconforming 
element for the second order elliptic eigenvalue problem. Let II be defined as in Section [O 
(up to some corresponding modification on the boundary for the boundary condition). We 
define 

(D 1) A a(Jlu h ,Ilu h ) 

h {pTlu h , Ilu h ) L 2 {n) ' 

Let A be the exact eigenvalue of the corresponding eigenvalue problem. We have 
(D.2) \ h > A. 
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This is to say A^ is an upper bound of A. Next, we shall show that A\ — A has the same 
convergence rate as A ft — A. In fact, we have [El Lemma 4.1], 

mq v f ^ _ a ( u - Uu h, u - u Uh ) \\p 1/2 (u - nn fe ) ||| 2(n) 

( j h ll^n^ll 2 ^ llP^nu*!!^ " 



Then estimates in (1C6|) and (1C23|) state that 



(D.4) 1 1 Ilit — Uh\\h % ||w — Uft || h + hi§ n or der term. 

Therefore we have 

(D.5) X h - A < /i 2s 

provided that it G H m+S (n) with < s < 1. 

With both the lower and upper bounds of the eigenvalue, we can combine them to get a 
more accuracy approximation, see |3H] for more details. 

One may wonder whether we can first compute the approximate eigenvalue A^ and its 
corresponding eigenf unction, say, w^, by the conforming finite element method. Then we 
use the canonical interpolation n nc u^ of some nonconforming space to define 



(D.6) A 



| p l/2 n n Cw 



h\\L 2 (n) 



Can we prove that A^ c is a lower bound of A? In general, the answer is negative unless 
the conforming method is of higher order than the nonconforming method, which is con- 
firmed by both the numerical and theoretical results in [4TJ]. This shows the superiority of 
nonconforming methods over their conforming counterparts. 
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